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Abstract 

We revisit a two right-handed neutrino model with two texture zeros, namely an indirect 
model based on A^ with the recently proposed new type of constrained sequential domi- 
nance (CSD2). involving vacuum alignments along the (0, 1, — 1) T and (1, 0, 2) T directions 
in flavour space, which are proportional to the neutrino Dirac mass matrix columns. In 
this paper we construct a renormalizable and unified indirect A^ x SU(5) model along 
these lines and show that, with spontaneous CP violation and a suitable vacuum alignment 
of the phases, the charged lepton corrections lead to a reactor angle in good agreement 
with results from Daya Bay and RENO. The model predicts a right-angled unitarity tri- 
angle in the quark sector and a Dirac CP violating oscillation phase in the lepton sector 
of S 130°, while providing a good fit to all quark and lepton masses and mixing angles. 
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1 Introduction 



The lepton mixing angles have the distinctive feature that the atmospheric angle #23 and the 
solar angle #12, are both rather large PQ. Direct evidence for the reactor angle #13 was first 
provided by T2K, MINOS and Double Chooz [2H1|. Subsequently Daya Bay [5], RENO [BJ, 
and Double Chooz [7] Collaborations have measured sin 2 (2^13): 

Daya Bay: sin 2 (26>i 3 ) = 0.089 ± 0.011(stat.) ± 0.005(syst.) , 
RENO: sin 2 (2(9i3) = 0.113 ± 0.013(stat.) ± 0.019(syst.) , (1.1) 

Double Chooz: sm 2 (2e 13 ) = 0.109 ± 0.030(stat.) ± 0.025(syst.) . 

This rules out the hypothesis of exact tri-bimaximal (TB) mixing [8], and many alternative 
proposals have recently been put forward [9], although there are relatively few examples which 
also include unification |10H13j . For example, an attractive scheme based on trimaximal 
(TM) mixing remains viable |14j , sometimes referred to as TM2 mixing since it maintains the 
second column of the TB mixing matrix and hence preserves the solar mixing angle prediction 
sin 612 ~ 1 /\/3. However there is another variation of TM mixing which also preserves this 
good solar mixing angle prediction by maintaining the first column of the TB matrix, namely 
TMi mixing [15]. 

Although there were models of TM2 mixing which can account for the smallness of the 
reactor angle [16], the first model in the literature for TMi mixing, which also fixed the 
value of the reactor angle, was proposed in [T7]. The model discussed in [T7] was actually 
representative of a general strategy for obtaining TMi mixing using sequential dominance 
(SD) [18] and vacuum alignment. The strategy of combining SD with vacuum alignment is 
familiar from the constrained sequential dominance (CSD) approach to TB mixing [19] where a 
neutrino mass hierarchy is assumed and the dominant and subdominant flavons responsible for 
the atmospheric and solar neutrino masses are aligned in the directions of the third and second 
columns of the TB mixing matrix, namely (4>\) oc (0, 1, — 1) T and (cj)^} oc (1, 1, 1) T . The new 
idea was to maintain the usual vacuum alignment for the dominant flavon, {(f>\) oc (0, 1, — 1) T 
as in CSD, but to replace the effect of the subdominant flavon vacuum alignment by a different 
one, namely either ((^120) oc (1,2,0) T or (^102) oc (1,0, 2) T , where such alignments may be 
naturally achieved from the standard ones using orthogonality arguments. 

We referred to this new approach as CSD2QI and showed that it leads to TMx mixing and 
a reactor angle which, at leading order, is predicted to be proportional to the ratio of the solar 
to the atmospheric neutrino masses, #13 = ^- — §■. The model was proposed before the results 
from Daya Bay and RENO, and the prediction turned out to be rather too small compared 
to the results in Eq. (jl.ip . More generally it has been shown that any type I see-saw model 
with two right-handed neutrinos and two texture zeros in the neutrino Yukawa matrix (as in 
Occam's razor) is not compatible with the experimental data for the case of a normal neutrino 
mass hierarchy [21]. However this conclusion ignores the effect of charged lepton corrections, 
and so an "Occam's razor" model which includes such corrections may become viable. 

In the present paper we construct a fully renormalisable unified x SU(5) model in 
which the neutrino sector satisfies the CSD2 conditions, and show that, with spontaneous CP 

It is interesting to compare the predictions of CSD2 to another alternative to CSD that has been proposed 
to account for a reactor angle called partially constrained sequential dominance (PCSD) [20]. PCSD involves 
a vacuum misalignment of the dominant flavon alignment to (e, 1, — 1) T , with a subdominant flavon alignment 
(1, 1, 1) T , leading to tri-bimaximal-reactor (TBR) mixing [2D] in which only the reactor angle is switched on, 
while the atmospheric and solar angles retain their TB values. 
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violation and a suitable vacuum alignment of the phases, the charged lepton corrections can 
correct the reactor angle, bringing it into agreement with results from Daya Bay and RENO. 
We shall use here similar techniques as in [22], where spontaneous CP violation with flavon 
phases determined by the vacuum alignment was discussed for the first time, in order to ensure 
that the charged lepton mixing angle correction (typically about ~ 3°) adds constructively 
to the 6i 3 angle from the neutrino sector (typically about ~ 5° — 6°) leading to #13 ~ 8° — 9°, 
within the range of the measured value from Daya Bay and RENO. In fact the present model 
is more ambitious, since it describes all quark and lepton masses and mixing angles, including 
predictions for all the CP violating phases. 

We demonstrate the viability of the model by performing a global fit to the charged 
lepton masses and the quark masses and mixing parameters. For the neutrino mixing angles 
we make a parameter scan and find very good agreement with the experimental data. We 
emphasise that the present A4 x SU (5) model represents one of the first unified "indirect" 
family symmetry models in the literature that has been constructed to date that is consistent 
with all experimental data on quark and lepton mass and mixing parameters where "indirect" 
simply means that the family symmetry is completely broken by the vacuum alignment H 

We emphasise that the idea of spontaneous CP violation has a long history |23j . However, 
in explicit flavour models using this idea only the positions of the phases in the mass matrices 
was predicted, but not the phases of the flavon fields themselves (see e.g. (23]). Spontaneous 
CP violation with calculable flavon phases from vacuum alignment was first discussed in [22] 
and demonstrated in example models based on A4 and S4. In this paper we shall use a 
similar approach where the A4 model is formulated in the real SO (3) basis (see e.g. [25] ) 
and where we only consider the real representations 1 and 3. In such a framework, one 
can either use a "simple" CP symmetry under which the components of the scalar fields 
transform trivially as <j>i — )■ <f>*, or a "generalised" CP symmetry which intertwines CP with 
A4 (see e.g. [25] and references therein). In the latter case, in our basis, the triplet fields 
would transform as 0, — > Us4>*, where U3 interchanges the second and third component. 
When complex 1' and 1" representations are used in a model, the U3 transformation then 
takes care of the fact that under CP the two complex singlets are interchanged with each 
other. However, as already mentioned above, in our model this will make no difference. CP 
symmetry leads to real coupling constants in a suitable field basis (after "unphysical" phases 
have been absorbed by field redefinitions). CP is subsequently spontaneously broken by the 
flavon vacuum alignment, which is controlled by additional Abelian symmetries Z3 and Z4, 
resulting in calculable complex flavon phases as in [22]. 

The layout of the rest of the paper is as follows: in the next section we discuss the general 
strategy we will adopt in our model. After a brief review of CSD2 we discuss charged lepton 
sector corrections to TMi mixing before we describe the method which we use to fix the 
flavon vevs. In section [3] we describe our model, the field content and symmetries and the 
resulting Yukawa and mass matrices. The justification for the chosen vacuum alignment 
including phases is given in section [3J In the subsequent sections we comment on the Higgs 
mass and then we give the numerical results from our global fit and scans. In section [7] we 
summarize and conclude and in the appendix we define our notations and conventions and 
give the messenger sector of our model. 

2 In fact the only other example of a unified indirect model with a realistic reactor angle that we are aware 
of is the last paper in [T5] based on Pati-Salam unification, however that model predicts an atmospheric angle 
in second octant. 
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2 The strategy 



Let us now describe our general idea in somewhat more detail, before we present an explicit 
GUT model example in the next section. As outlined in the introduction, we are combining 
three ingredients which finally result in a highly predictive unified flavour model. These 
ingredients are: 

• CSD2 for the neutrino mixing angles 0^-, 

• charged lepton mixing contributions as they are typical in GUTs, 

• spontaneous CP violation with aligned phases. 

We now briefly describe these three concepts and the resulting new class of models. 



CSD2 in the neutrino sector 

In models with CSD2 [T7], the neutrino mass matrix is dominated by two right-handed 
neutrinos with mass matrix Mr = diag(M J 4, M#) and couplings to the lepton doublets 
A = (0, a, —a) T and B = (6, 0, 2b) T d such that the neutrino Yukawa matrix takes the form 
Y v = (A,B), in left-right convention. A summary of the used conventions is given in the 
Appendix. 

After the see-saw mechanism is implemented, CSD2 leads to the following light effective 
neutrino Majorana mass matrix: 



M u = m a 1 -1 \+m b = m a 1 -1 , (2.1 
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where m a = ^ftp-, m b = ~m~b~' an< ^ wnere a * s ^ ne relative phase difference between m a and m&. 
We define e = |m;,|/|m a |, and assume e <C 1 leading to a normal mass hierarchy in accordance 
with SD. As discussed in Appendix [A] we use here different conventions than in the original 
CSD2 paper [T7] which are more convenient in the context of SU(5) GUTs. 

Only three parameters, e.g. m a , e and a, govern the neutrino masses and mixing param- 
eters. For the mixing parameters, the predicted values are, to leading order in e (from |17] 
with adapted conventional): 

1 e 5 
^ 3 ~ "71 ~ \7I C ° S ' <5i3 « vr - a + e-sina , (2.2) 

, «2 ~ —ol + 2e sin a , (2-3) 

V2 

The mixing scheme resulting from CSD2 can be identified as trimaximal mixing of type 1 (i.e. 
TMi [H]) but with a predicted value of the neutrino 1-3 mixing, 0^ 3 = ~ 5° - 6°. With 

neutrino mass m\ = 0, only one Majorana CP phase is physical. Without charged lepton 
corrections, would be identical with the leptonic Dirac CP phase 5. Let us also note that 
CSD2 predicts a deviation of from 45°, depending on the phase a. 

:i B = (b, 2b, 0) T was also considered in [T7j, but here we shall not consider it further. 
4 Compared to the notation of [17], we have changed, for instance, a — > —a. 
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Charged lepton mixing contribution in GUTs 

In GUT models the charged lepton Yukawa matrix is generically non-diagonal in the flavour 
basis, due to the close link between the charged lepton and the down- type quark Yukawa 
matrices, which typically provides the main origin of the flavour mixing in the quark sector. 
With the Cabibbo angle 9c being the largest mixing in the quark sector, the mixing in Y e is 
often dominated by a 1-2 mixing 9\ 2 as well, such that the relevant part of (the hierarchical 
matrix) Y e can be written as 

/0 ce i/3 0\ 

Y e « * d , (2.5) 
\0 * */ 

where c, d, and /3 are real and where the entries marked by a V are not relevant for our 
discussion here. With c <C d one can read off to leading order the values for the complex 1-2 
mixing angle (for more details see also Appendix [AJ that are 



I'll 



12 \-(3 + tt forc/d<0 



Since we will have c/d < in our example GUT model in the next section, let us consider 
this case also in the following discussion. 

In explicit GUT models, 9f 2 is typically related to the Cabibbo angle by group theoretical 
Clebsch factors from GUT breaking, as has been discussed recently, e.g. in jl(J|,lll j . In many 
GUT models, in particular in those where the muon and the strange quark mass at the GUT 
scale is predicted by such a Clebsch factor as m^/mg = 3 [27,, but also if the Yukawa matrices 
Y e (and Yd) are (nearly) symmetric with a zero in the (0,0)-element [10] . 9\ 2 is predicted as 



9r 

??2 « ■ (2-7) 



3 

In the example GUT model in the next section we will see explicitly how such a prediction 
arises in an SU(5) GUT. 

The leptonic mixing parameters, defined via C/pmns = U e ut, are a combination of the 
mixing from the neutrino and the charged lepton sectors. Making use of the fact that, to 
leading order, 9 23 = 9f 3 = 5\ 2 = 5 23 = 0, and using the CSD2 expressions from above for the 
neutrino sector, and general formulae for the charged lepton mixing contributions of j!9]l28]l29j 

(2.8) 
(2.9) 
(2.10) 



S23e" 


i #23 


— s 23 e 


i% - 




si3e" 


i5i3 


= 0? 3 e _ 


i% - 




si2e~ 


i<5i2 


= si 2 e~ 


i% _ 





we obtain (up to 0(e)) 



? 23 pa 45° - e cos a , (2.11) 
? 13 «-^-cos(/3-a)^|, (2.12) 

ne 

?12 ~ 35.3° + cos p , (2.13) 
v2 
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where e ~ ~ 8.4°. When the phases a and (3 are fixed by the vacuum alignment, and 
when also 9f 2 is predicted from the GUT structure, as both will be the case in our model, 
all three mixing angles and also the CP phases 5 and a2, are predicted. Thus, the resulting 
models of this type can be highly predictive. 

We would like to note here already that in the explicit GUT model in the next section, 
we will construct a vacuum alignment such that a = vr/3, leading tcH 

#23 ~ 45° - I « 41° , (2.14) 

close to the best fit value for the normal hierarchy case from global fits to the neutrino 
data [30]. The alignment of /3 will satisfy /3 = a + tt, such that the neutrino and charged 
lepton contributions to #13 simply add up, leading to (with 9f 2 = #c/3) 

in agreement with the recent measurements. With these values of a and (3, it also turns out 
that 612 is predicted somewhat smaller than 35°, namely 



'12 



33° . (2.16) 



This value of 612 could be distinguished from the tribimaximal value by a future reactor 
experiment with ~ 60 km baseline |31j . 



Spontaneous CP violation with aligned phases 

Finally, the third ingredient is spontaneous CP violation with aligned phases of the flavon 
vevs, using the method proposed in [22]. To give a brief summary of this method, let us note 
that phase alignment can very simply be achieved using discrete symmetries when the flavon 
vevs effectively depend on one parameter, i.e. when the direction of the vevs is given by the 
form of the potential. This remains true even in the presence of "generalised" CP transforma- 
tions as long as these CP transformations fix the phases of the involved coupling constants. 
Working example models with A4 and S4 family symmetry can be found in |22j . Note that 
S4 is in agreement only with "simple" CP, while the "generalised" CP transformation for A4 
interchanges the complex singlet representations [26] . In both cases all the coupling constants 
are forced to be real in a suitable field basis. 

To illustrate the phase alignment, let us consider a case with a flavon field £ which is a 
singlet under the family symmetry and singly charged under a Z n shaping symmetry (with 
n > 2). Then typical terms in the flavon superpotential, which "drive" the flavon vev non- 
zero, have the form 

P (7^2 TM^. (2.17) 

The field P is the so-called "driving superfield", meaning that the F-term |.Fp| 2 generates the 
potential for £ which enforces a non-zero vev. A is the (real and positive) suppression scale 

5 We note that the choice a = tt/3 is motivated by the current data which favours #23 in the first octant. 
On the other hand, one can in principle also construct other models with different values of a, and there are 
also other options for j3 and 9\ 2 , which may lead to interesting alternative models. In this sense, the strategy 
described here leads to a whole new class of possible models. 



5 





Ti 


T 2 




F 


N x 


N 2 


#5 


H 5 


H45 


H45 


H24 


S 


nrr/r\ 

bU(o) 


10 


10 


10 


_ 
5 


1 


1 


5 


5 


45 


45 


24 


1 


Aa 

1 


1 


1 


1 


3 


1 


1 


1 


1 


1 


1 


1 


1 


U(1)r 


1 


1 


1 


1 


1 


1 

















2 


Z 4 


3 


3 


3 








2 


2 








2 


1 


2 


z 4 


3 


3 


3 





2 


2 


2 





2 





1 


2 


z 3 


1 


2 








1 


2 














2 





z 3 


1 


1 








2 











1 


2 


2 





z 3 





2 


2 


1 





2 


2 





1 


1 


2 


1 


z 3 











2 














1 


2 


1 






Table 1: The matter and Higgs fields in our model and their quantum numbers. 



of the effective operator, and M here is simply a (real) mass scale. From the potential for £, 

2 



s =fM 2 



the vev of £ has to satisfy 



e 



A n-2 



± A n-2 M 2 



Since the right side of the equation is real, we obtain that 



arg«£» 




. ,n 
l,...,n 



for "-" in Eq. (j2T7|) . 
for "+" in Eq. ([2~TTjl . 



(2.18) 



(2.19) 



(2.20) 



For example, with a Z 3 shaping symmetry and a "+" in Eq. (|2,17p . only multiples of 7r/3 
are allowed for arg ((£)). We will use this method for the relevant flavons to constrain their 
phases. In the ground state, one of the vacua (with a fixed phase) is selected, which finally 
determines also the two phases a and /3 relevant for the predictions in the lepton sector. 

Furthermore, we note that we will also use the phase alignment to generate the CP 
violation in the quark sector, predicting a right-angled unitarity triangle, which is in excellent 
agreement with the present data (making use of the quark phase sum rule from [32]). 

We now turn to an explicit GUT model, where the above described strategy is applied. 



3 The model 

In the following we will construct an Aa x SU (5) model with CSD2 |17| in the neutrino sector. 
The model follows the strategy described in the previous section, such that the charged lepton 
mixing contribution to #i 3 adds up constructively with the 1-3 mixing in the neutrino sector 
to 612, ~ 8° — 9°, with the phases fixed by the "discrete vacuum alignment" mechanism |22j . 

The matter and the Higgs sector of the model is summarised in Table [1] while the required 
flavons are shown in Table [2j The superpotential after integrating out the heavy messenger 
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SU{5) A 4 U(l) R Z 4 Z 4 Z 3 Z 3 Z 3 Z 3 
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Table 2: The flavon field content of our model. 



fields, see Appendix [Bl and suppressing order one coefficients reads 

nv = 6^i 2 + kivf , (3.i) 

W v = ]-( J ff5i ? )(023iVi) + j {H 5 F){cj> W2 N 2 ) , (3.2) 

m = ^3^2^5^(^102)^24 + ^e i02 H 5 F(T 2 ^ W2 )H 24 + ^F(T 2 02 3 )^45^24 + jH 5 F(T 3 (l> 3 ) , 

(3.3) 

W u = plf ff 5 ^i + ^T x T 2 H^l + -^TiH^l + jT 2 T 3 H 5 h + T 3 2 ^ 5 , (3.4) 

where A denotes the messenger scale. The flavon potential, which gives rise to the vevs of the 
fields 0j, £j and 9{ will be discussed separately in the next section. Note that the flavons of 
type which enter the Yukawa couplings will be aligned with real vevs while the flavons of 
type 9 and £ will generally acquire complex vevs with precisely determined phases. The above 
superpotential gives rise to the flavour structures in the neutrino sector, in the down-type 
quark and charged lepton sectors, and in the up-type quark sector. 

Neutrino sector: From the flavon potential, to be discussed in the next section, the two 
triplet flavons entering the neutrino Yukawa sector are aligned along the directions 




23/ ~ 1 , Wl02j 



(3.5) 



where both alignments are real. Inserting the above vacuum alignments, the real vev and 
the vev (£2) with a phase of —tt/3 into the superpotential leads to a Dirac Yukawa matrix 
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and a right-handed heavy Majorana mass matrix of the form: 



b 
a 
-a 2b, 



where Ma, a and b are real and Mb has a complex phase of —tt/3. The (type-I) seesaw 
formula leads to a simple effective light neutrino mass matrix of the form given in eq. (|2,lj) 
where the relative phase difference a between m a and nif, is now fixed to be ir/3. This form 
of M v gives #^3 ~ 5° — 6° for the 1-3 mixing in the neutrino sector, which will finally add up 
with the charged lepton mixing contribution. 

Down-type quark and charged lepton sector: Turning to the down quark and 
charged lepton sector, two further triplet flavons enter: 



(3.7) 



where (^2) is aligned to be real. The phase of ((^3) will turn out to be unphysical. Furthermore 
the singlet 62 with a phase of ir/2 and the singlet #102 with a phase of 4ir/3 enters. Plugging 
in the vevs of the flavon fields leads to the following structure of the Yukawa matrices (in 
left-right convention) for the down- type quarks and charged leptons: 

ie 2 \ / -3/2weio2 

Y d =\Oe W2 e 23 2weio2-e 2 3 and Y e = -3/2 i e 2 9/2 e 23 0| , (3.8) 
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where ui = e i7T 1 / 3 , cf . section [2j The e, are proportional to the order one couplings which we 
have not written down explicitly and possible Higgs mixing angles. For the sake of simplicity 
we only show here the proportionality to the dimensionful quantities 

£2 ~ ^|(6 , 2)(02)|, £102 ~ ^ I (0102) (0102) I, ^23 ~ ^yKfe)!, ^3 ~ T^l(03)|, (3.9) 

where ^24 is the vev of H24 ■ We also note that we do not use the common Georgi-Jarlskog 
relation m^/m s = 3 [27] at the GUT scale but rather m^/rris = 9/2 [33102]. The reason for 
this is that recent lattice results, see, e.g. [35] suggest a much smaller error for the strange 
quark mass than the PDG quotes. And since we are in the small tan ft regime and no large 
SUSY threshold corrections can correct the second generation GUT scale Yukawa coupling 
ratios we have to use the more realistic relation mentioned above. Explicitly, from the vevs 
of #24 and H$ we get a relative factor of —3/2 for e 2 and €102 and the 9/2 from H24 and 
For the third generation we use b — r Yukawa unification which is possible for small tan /3 due 
to the large RGE effects induced by the top mass. 

For the 1-2 mixing in the charged lepton sector, we nevertheless obain 9f 2 ~ 0c/3, where 
9q ~ 0.23 is the Cabibbo angle. The corresponding phase 5\ 2 is chosen (see section [2] and 
appendix |A] for conventions), such that the charged lepton mixing angle correction 9\ 2 is in 
phase with the neutrino reactor angle 9^ and the two angles add together constructively to 
yield the physical reactor angle #13. 
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Up-type quark sector: Finally the up-type quark sector only involves singlet flavons 
with real vevs and gives a real symmetric Yukawa matrix of the form, 





K 





K 











t ' u 



Y u = \b u c u d u , (3.10) 



where the dependence on A and the flavon vevs reads 

11 A 1 ' u W ' " A2~' 14 A - ' ^ ' 

Note that e u is coming from a renormalisable coupling. The zero texture in the quark sector 
means that we can successfully apply the quark phase sum rule of [32] due to our choice of 
phases. 



4 The vacuum alignment 

We have in total seven flavon fields which transform as triplets under A4, see Table [2J pointing 
in the following directions in flavour space, 



(4.1) 




211/ ~ 1 , Will) ~ 1 , \CP231 




(4.2) 



Apart from (</>i) and ((^3), the vevs of the above listed flavons will be aligned real using the 
phase alignment mechanism proposed in [22]. The phases of {(pi) and (^3) have no physical 
implications and hence will be set real for definiteness. The first three vevs form a basis in 
flavour space, while the second three alignments are proportional to the (real) columns of the 
tri-bimaximal mixing matrix. In our model, instead of 0m (which is used in the CSD [18^19] 
models), we require the following (real) alignment, 



>102 



(4.3) 



in the neutrino sector, similar to a recently proposed flavon alignment [T7] but with the phase 
fixed as explicitly shown and discussed below. 

The principal assumption of our model is that CP is conserved above the flavour breaking 
scale, and is spontaneously broken by the CP violating phases of flavon fields. With this 
assumption we can not only reproduce the correct mixing angles but can also make definite 
testable predictions for the CP violating phases in the lepton sector. In order to do this we 
will fix the phases of the following flavon vevs to 

am = , a 2 n = , a 23 = , a 2 = , a W 2 = , (4.4) 
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1 




1 


3 


2 


2 


2 
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A 2 


1 


3 


2 





2 














A 3 


1 


3 


2 


2 


2 











1 


Am 


1 


3 


2 


2 


2 














P 


1 


1 


2 





















Table 3: The driving field content of our model. Note that we only show here one P field. 
Indeed one has to introduce as many P fields as operators to fix the phases of the flavon 
fields. Since they will have all the same quantum numbers they will mix and we can go to a 
basis where the terms to fix the phase for each flavon is separated from the others. This was 
discussed in the appendix of |22j . 



where a, stands for the phase of (</>«). Furthermore we have some singlet flavons with non- 
vanishing vevs of which some will have non-trivial phases. In this choice we have also ignored 
possible signs which means that the phases are fixed up to ±7r. We can fix the phases by 
using appropriate Z n shaping symmetries as described in our previous paper [22], see also 
section [2j 

The method can be understood easily for the singlet flavon vevs. Their superpotential 
reads 

W = ^(£ 3 - M 3 ) + + M 3 ) + j(f u - M 3 ) 

+ ^(91 - M 4 ) + ^-(9f 02 - M 3 ) + ^(p 3 02 - M 3 ) + ^ - Af 3 ) , (4.5) 

where M is a generic mass scale which we assume to be positive. The list of the driving fields 
is given in Table The F terms for P will then fix the flavon vevs of the singlets up to 
a discrete choice. Note that for the sake of simplicity we have only introduced one P field. 
Indeed, we need one P field for every singlet. Since they all have the same quantum numbers 
they will mix and we can go to a basis where all terms are disentangled as in the equation 
above, see the appendix of [22]. For the singlet flavons here we choose and (^102,23) to 

be real, (#2) to be imaginary, (#102) to have a phase of 4ir/3 and (£ 2 ) to have a phase of — tt/3. 

We come now back to the phases of the triplet flavon vevs which can be fixed in the same 
way after the direction in flavour space is fixed. Note that the phases ot\ and a 3 are not 
fixed in our model. This is also not necessary. The flavon (j>i does not couple to the matter 
sector and hence its phase does not appear in the mass matrices. It will only be used in 
orthogonality relations where the phase of the vev does not matter. The flavon (^3 couples 
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nevertheless to the matter sector. But as we have seen before it determines the 3-3 element 
of the down-type quark and charged lepton Yukawa matrix and its phase can be absorbed in 
the right-handed fields such that this phase renders unphysical. 

In this section we will use an explicit notation for the contraction of the A4 indices. We 
use the standard "50(3) basis" for which the singlet of 3 (g) 3 is given by the S"0(3)-type 
inner product The two triplets of 3 <g> 3 are constructed from the usual (antisymmetric) 
cross product 'x' and the symmetric star product V (see, for example, [25]). 

We start with the alignment of the triplet flavons 0j, i = 1, 2, 3, which can be aligned via 

W = Ai ■ (fa * <fH) + O i;i (& • 0,) + ^ ((0 2 • 2 ) 2 - M 4 ) . (4.6) 

Solving the F-term conditions of Ai aligns the flavons in one of the three standard directions 
and the F-term conditions of Oi-j makes them orthogonal to each other. By convention we 
let them point in the directions as given in eq. (|4.ip . For «2 we choose the value {a\ and 
03 remain undetermined). In Appendix [B] we will discuss the messenger sector of our model. 
After integrating out heavy messenger fields we end up only with the effective operators 
written here and in the following. 

We now turn to the flavons 023 , 0m arid 02ii : F° r 0m we use a slight modification of 
the alignment in the recent SU(5) x T' model [T3] without auxiliary flavons, 

p 

W = A ni ■ (0ni * 0ni + 0111P1H + 0mpiii) + --2 ((0iii ' 0m) 2 - M 4 ) 



A 2 

+ ^ (/'in + PniPni + Pin ~ M ") ■ ( 4 - 7 ) 



It gives the desired alignment and (0m) can be chosen to be real. 
Starting from this the other two alignments can be realised by 

W = Ol;23(01 • 023) + Oin ; 23(0111 ■ 023) + Olll;211 (0111 ■ 021l) + O 2 3;21l(023 " 021l) 

P P (4.8) 

+ -r ((0211 * 021l) • 0211 - M 3 ) + - ((023 • 023)P23 ~ M 3 ) . 

The orthogonality gives the desired directions and (02ii) can be chosen to be real. The phase 
of (023) is a bit peculiar. Above we have fixed (P23) to be real and hence also (023) can be 
chosen to be real. In the first operator the vev of 0i enters again and independent of the 
phases a (0, 1, —1) alignment is always orthogonal to a (1, 0, 0) alignment. 

Now we have everything together for the last missing non-trivial alignment 

p 

VV = O211;1O2(01O2 • 021l) + O2;1O2(01O2 • 02) + "T" ((0102 " 0102)^102 ~ M 3 ) . (4.9) 

The direction is again fixed by orthogonality conditions. The vev of 0io2 can be chosen to be 
real (remember that also (/O102) is real). 



5 The Higgs mass 

In our model we assume b — t Yukawa coupling unification at the GUT scale. This happens 
in the MSSM only for large tan (3 via SUSY threshold corrections or small tan f3 due to large 
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;T/dof = 2.05/3 



+2(T 
+ \<T 

Oo- 

-\(T 



T 



-la 



1 "12 "13 "23 °CKM 



Figure 1: Pictorial representation of the deviation of our predictions from low energy ex- 
perimental data for the charged lepton Yukawa couplings and quark Yukawa couplings and 
mixing parameters. The deviations of the charged lepton masses arc given in 1% while all 
other deviations are given in units of standard deviations a. 



RGE corrections by the top mass. We have decided for the second solution such that we can 
also neglect SUSY threshold corrections in our fit later on. 

Nevertheless, the MSSM with small tan /3 prefers very light Higgs masses which is in 
conflict with the recent discovery of a Higgs-like particle with a mass of about 126 GeV [36] , 

A possible solution to this problem is given by the NMSSM, for a review see [37j where the 
Higgs can have the right mass even for small tan/3. In fact our symmetries forbid a /i-term 
because the combinations and are charged under the shaping symmetries. But 

we have checked that we can add a singlet field S which couples simultaneously to this two 
combinations. For convenience we have listed the field S in Table [TJ 

An explicit S 3 term in the superpotential is forbidden in the limit of unbroken U(1)r sym- 
metry (i.e. before SUSY breaking) and by the shaping symmetries but is needed to stabilize 
the Higgs potential in the scale invariant NMSSM. But we note that there are still various 
possibilities to stabilize the potential for S. This could be done, for instance, by introduc- 
ing an additional U(l)' gauge group where the potential is stabilized by the U{1)' D-terms. 
For a description of this and references, see the review article [37] • We only note that it is 
straightforward to introduce auch a U(l)' in our model by charging the Higgs and matter 
fields appropriately which does not alter the flavour sector. Alternatively, the 5 3 term could 
be generated non-perturbatively, breaking the shaping symmetries in an F-theory framework, 
see, for instance, [38]. We will not go here into more detail on this model building aspect and 
only like to note that our flavour model is compatible with some NMSSM variants and hence 
we can have a realistic Higgs mass. 

6 The fit and numerical results 

Here we will present the results of a numerical x 2_n t °f the high energy parameters of the 
Yukawa matrices to the low energy charged lepton and quark masses and quark mixing pa- 
rameters. Afterwards we will present the predictions for neutrino masses and mixing. 

For the RGE running of the Yukawa matrices we have used the REAP package [39J and 
calculated with it the masses and mixing angles at low energies. Note that we have used 
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P PI V A TY1 PT PT" 


\f a 1 ii p 


a u 


— 3.01 ■ 10 


K 


-2.66 • 10" 4 


( 'u 


-2.57- 10~ 3 


d n 


3.09 • 10~ 2 


6 u 


2.05 




-3.57 • 10" 5 


£102 


3.17- 10~ 5 




1.62 • 10~ 4 




1.24 • 10" 2 


tan p 


1.49 



Table 4: Values of the effective parameters of the quark and charged lepton Yukawa matrices 
and tan/3 for Msusy = 750 GeV. The numerical values are determined from a x 2_n t to 
experimental data with a x 2 per degree of freedom of 2.05/3. 



Quantity (at mt(mt)) 


Experiment 


Model 


Deviation 


y T in 10~ 2 
y^ in 10~ 4 
Ve in 10~ 6 


1.00 
5.89 
2.79 


1.00 
5.89 
2.79 


-0.277 
0.097 
-0.016 


Vb in 10~ 2 
y s in 10~ 4 
Vd in 10~ 6 


1.58 ±0.05 

2.99 ±0.86 
K q+6.8 


I. 64 
2.95 

II. 7 


1.088 
-0.226 
-0.639 


yt 

y c in 10~ 3 
Vu in 10~ 6 


0.936 ±0.016 
3.39 ± 0.46 
7 01+ 2 - 76 


0.939 
3.40 
7.59 


0.159 
0.223 
0.209 


flCKM 

ziCKM 
"23 
nCKM 
a 13 

<5cKM 


O9K7+0.0009 

u.zzo(_ 0010 

041 ^+ a0011 

U.U41O_ 0012 

0.0036 ± 0.0002 

1 9n9Q+ - 0786 

i.zuzo_ 00431 


0.2257 

0.0409 
0.0036 
1.1975 


0.026 

-0.488 
-0.002 
-0.113 



Table 5: Fit results for the quark Yukawa couplings and mixing and the charged lepton 
Yukawa couplings at low energy compared to experimental data. The values for the Yukawa 
couplings are extracted from [5T] and the CKM parameters from [55]. Note that the experi- 
mental uncertainty on the charged lepton Yukawa couplings are negligible small and we have 
assumed a relative uncertainty of 1 % for them. The \ 2 P er degree of freedom is 2.05/3. A 
pictorial representation of the agreement between our predictions and experiment can be found 
as well in Fig. [TJ 
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the RGEs of the MSSM. Possible RGE effects due to including a variant of the NMSSM are 
neglected. On the one hand we can expect this effect to be flavour blind leading only to a 
rescaling of the GUT scale parameters and on the other hand, in the scale-invariant NMSSM 
for example, the RGE effects come from the coupling A which can be expected to be small 
from recent global fit results [30]. Furthermore, SUSY threshold corrections are negligibly 
small due to the small tan (3 and hence are not included in the fit. 

For the charged lepton and quark masses and their errors at the top scale mt(mt) we 
have taken the values from [3T] and for the CKM parameters the PDG values [12]. Note 
that the experimental errors for the charged lepton masses are tiny and we have estimated 
the theoretical uncertainty from higher order effects to 1 %, and we will assume this as their 
errors instead. 

The Yukawa matrices depend on nine real parameters (five from the up-type quarks and 
four from the down- type quarks and charged leptons) . Furthermore we have included tan /3 
as a free parameter in the fit. The unification of the b and the r Yukawa coupling at the 
GUT scale depends strongly on this parameter. On the contrary, the masses and mixing 
angles depend only very weakly on the SUSY scale which we have therefore fixed to Msijsy = 
750 GeV. 

The fit results are summarised in Figure [1] and Tables 0J and [5j We have fitted ten 
parameters to thirteen observables with a y 2 °f 2.05 and hence we can say that our model 
describes the data very well§ 

We turn our discussion now to the neutrino sector. Here we did not fit the parameters to 
the observables because here we are more interested in the allowed ranges and correlations 
between different observables which help in distinguishing this model from other models. 

The effective neutrino mass matrix from eq. (j'i.ip depends on three parameters. The 
neutrino mass scale m a the perturbation parameter e and the relative phase a. The phase 
a in our model is 7r/3 as discussed in section [2J Hence, only two real parameters m a and e 
completely determine all observables in the neutrino sector. 

We have varied these two parameters randomly and the results are shown in Figure [2] 
where we have used as constraint the fit results of the Bari group [30] . The blue dots agree 
with all experimental data within 3<r while the red dots agree even within la. The dashed 
lines in the plots label the corresponding allowed ranges of the observables on the axes. The 
la range of the leptonic Dirac phase 5 is shown in black because it is not measured directly 
and the fit results should be taken with a grain of salt. In the scan we also did not include it 
as a constraint. 

We are everywhere in good agreement with the experimental data and we find clear 
correlations. Especially, noteworthy is the value for #23 which lies around 38.5°. We also make 
precise predictions for the CP violating phases. One of the Majorana phases is unphysical 
because one neutrino remains massless. The Dirac CP phase has a value of 5 ~ 130° and the 
physical Majorana phase is Q2 ~ 315°. The Jarlskog determinant Jqp is around 0.025 and the 
effective neutrino mass for neutrinoless double beta decay m ee is of the order of 3 x 10 -3 eV, 
beyond the reach of current experiments. 

6 We note that while we get an excellent fit for the quark masses themselves, as given in the PDG review, 
there is some tension with QCD results which favour y s /y<i ~ 19 [43] . while our fit yields ys/]Jd = 25.3. We 
remark that this tension is the same that one also gets with the more conventional GJ relation instead of the 
Clebsch factors 9/2 and 3/2 used here, so it is not particular for our model. In our fit, we have not included 
Vs/yd as constraints, but we would like to note that future even more precise results on the quark masses, 
including lattice results, can provide powerful additional constraints on unified flavour models. 
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Figure 2: The correlations between 6*13 and the other two mixing angles and the two physical 
phases in PCSD2. The regions compatible with the la (3a) ranges of the mass squared 
differences and the mixing angles, taken from [3D], are depicted by the red (blue) points and 
delimited by dashed lines in corresponding colours. The la region for the Dirac CP phase is 
shown in black. 
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7 Summary 



We have constructed a unified A4 x SU (5) model featuring the new type of constrained 
sequential dominance CSD2 proposed recently in [T7]. The A4 x 577(5) model, with the 
CSD2 vacuum alignments (0,1, 1) T and (1,0, 2) T , provides an excellent fit to the present 
data on quark and lepton masses and mixings, including the measured value of the leptonic 
mixing angle #13 from Daya Bay and RENO, with testable predictions for the yet unknown 
parameters of the leptonic mixing matrix. 

The main idea of the present model is that, with a strong normal hierarchical spectrum 
(with m\ = by construction since there are only two right-handed neutrinos) the 1-3 angle 
in the neutrino sector, Q\%, is related to a ratio of neutrino masses by 9^ = ^zfy, leading to 
#13 ~ 5° — 6°. In addition, the reactor angle receives another contribution from mixing in the 
charged lepton sector. The charged lepton mixing induces a correction to #13 of ~ 3° which 
adds up constructively with 9\^ to give 

#13 ~ 8° - 9° , (7.1) 

within the range of the measured value from Daya Bay and RENO. The constructive addition 
of the neutrino and charged lepton mixing angles is achieved by assuming high energy CP 
invariance which is spontaneously broken by flavon fields whose phases are controlled using 
Abelian Z3 and Z4 symmetries as proposed in [22]. We emphasise that in our approach one can 
either use a "simple" CP symmetry, under which the components of the scalar fields transform 
trivially as ^ (^*, or a "generalised" CP symmetry (see e.g. [26] and references therein) 
where, in our basis, the triplet fields would transform as — > U^cp*, with U3 interchanging 
the second and third component of a triplet representation. 

The resulting unified flavour model is highly predictive, as described in section [6] since 
only two parameters determine the neutrino mass matrix, while the charged lepton corrections 
are fixed by the GUT framework: In particular, for the Dirac CP phase 8, for the one physical 
Majorana CP phase oli and for the atmospheric angle #23 we obtain the predictions 

5 » 130° , a 2 ~ 315° and 9 23 ~ 38.5° . (7.2) 

The predictions for 5 and #23 will be tested by the ongoing and future neutrino oscillation 
experiments. In addition, for 6>i 2 , we predict a value of 

#12 ~ 33° , (7.3) 

which is slightly smaller than the tribimaximal mixing value but may be tested by a future 
reactor experiment with ~ 60 km baseline, which could measure 9yi with much improved 
precision [31]. Furthermore, in the quark sector, we obtain a right-angled unitarity triangle 
(with a ~ 90°) from the same vacuum alignment techniques for the phases [22], realizing the 
phase sum rule of [32]. 

In summary, we have presented a highly predictive new unified model for fermion masses 
and mixing, which, in fact, represents the first unified indirect family symmetry model in the 
literature that has been constructed to date that is consistent with all experimental data on 
quark and lepton mass and mixing angles, and makes definite predictions for CP phases in 
both the quark and lepton sectors. 
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A Conventions and Notations 

In this section we want to summarize briefly our conventions and define some notation used 
throughout the main text. We will follow mainly the notation of The only difference is 
a sign in the Majorana phases. 

The Yukawa couplings follow the left-right convention 

£ Yuk = -Yij¥A H + H - c - > (A.l) 
and for the effective light neutrino mass matrix we use the convention 

C = —UiM^jL] + H.c. , (A.2) 

where L is the lepton doublet. 

In the quark sector we define the CKM matrix by 

^ckm = U U U\ = R 23 U 13 R 12 , (A.3) 

where U u (Ud) is a unitary matrix diagonalising Y u Yu {YdX\) and 

Q12 s 12 e~ iSl2 0\ 
U 12 = I -s 12 e iS ^ cia , (A.4) 
1/ 

and similar for U 23 and Ui 3 . We use c\ 2 and s\ 2 as abbreviations for cos #12 and sin Q\ 2 . The 
matrices R 2 % and R± 2 are C/23 and Ui 2 with the complex phases set to zero. In this case #13 
coincides with the CKM phase <5ckm- 
For the PMNS matrix we use 

C/pmns = U e U$ = R 23 U 13 R 12 diag(e- iai / 2 ,e- iQ2 / 2 , 1) , (A.5) 

where the neutrino mass matrix is diagonalized via 

U u M v Mlul = diag(mf , m|, m§) . (A.6) 

and Uv = ( n °te the Hermitian conjugation). This conventions imply a complex 

conjugation of the neutrino mass matrix M v compared to our previous CSD2 paper [T7] and 
also the sign of the Majorana phases here is different. 
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Table 6: The messenger field content of our model. Every line represents a messenger pair 
which receives a mass larger than the GUT scale and no cross terms are allowed. In the main 
text we labelled the messenger mass scale generically with A. 
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Figure 3: The supergraphs before integrating out the messengers for the fiavon sector (only 
diagrams arc shown which give non-renormalizable contributions). 

B The Renormalizable Superpotential 

In this appendix we discuss the full renormalizable superpotential including the messenger 
fields which after being integrated out give the effective operators as discussed before. 

We start with the superpotential bilinear in the fields which is in our case only the mass 
terms for the messengers 

Wr = + M r .Tifi + M Si EiEi + M n ^A + Mr^f; . (B.l) 

The full list of messenger fields is given in Table [6] where every line is a messenger pair which 
receives a mass larger than the GUT scale so that they can be integrated out to give the 
desired effective operators. To simplify the notation before we have introduced the messenger 
scale A as shorthand which is related to the individual messenger masses with order one 
coefficients. 

Note that in the superpotential bilinear in the fields no /x-term for the Higgs fields appears. 
This term is forbidden by symmetries and in combination with a NMSSM like mechanism helps 
to increase the Higgs mass to the experimentally determined value. A possible singlet field 
S with couplings S(H^H^ + H^H^) would not appear anywhere else in the superpotential 
with the symmetries and field content as specified in Tables [HOE! and [6j 

The next step in our discussion of the renormalizable superpotential is the fiavon sec- 
tor. The full potential for this sector reads (dropping for the sake of simplicity order one 
coefficients) 

Wflavon = Ol;20102 + Ol;30103 + C>2;30203 + Olll;211 0111 0211 + Olll;230111023 
+ O23;2110230211 + O2 ; 1O20201O2 + O211 ; 1O2021101O2 + Ol;2301023 

+ ^10101 + ^20202 + ^30303 + Mil (0111 + 011lPlll + 011lPlll) 

+ Pr 9 £ u + f 9 f u + PTj + f 8 0| + f 8 8 2 2 + PT 2 7 + f 7 (0? n + p 2 in + p? u ) 
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Figure 4: The supcrgraphs before integrating out the messengers for the down-type quark 
and charged lepton sector. 

+ P<p2nr 6 + 0| n f 6 + p&t 5 + £ 2 2 r 5 + P£iT A + e?f 4 + Pp 23 r 3 + (4? + 4»)f 3 
+ Ppi 02 r 2 + (4) 2 + p? 02 )f 2 + P9 W2 r l + 9 2 102 f 1 . (b.2) 

The first three lines of this superpotential have already been discussed in the flavon alignment 
section H] while the last four lines are needed to fix the phases of the various flavon vevs. For 
instance, the messenger pair Ti and Ti gives after integrating out the effective operator 
l/AP0f O2 where in this case A stands for Mr 1 multiplied by real order one couplings. This 
operator fixes the phase of (6*102) up to a discrete choice as discussed before. 

We will not list here all of the effective operators because they have already appeared in 
our superpotential for the flavon alignment and they can also be read off from the diagrams 
in Figure [3] after contracting the messenger propagators to points. 

For the renormalizable couplings including the matter and Higgs fields we find the renor- 
malizable superpotential (again dropping order one coefficients) 

= T 3 i? 5 £ 3 + F0 3 S 3 + T 2 ^ 23 T 1 + i^TiHi + FH M B X + T 2 </> 102 T 2 + H 5 T 2 E 2 



+ #io 2 ^ 2 ^i + 7i(/> 2 T 3 + H 5 T 3 E 3 + # 2 H 3 Hi , (B-3) 

C = T i^ 3 + &n2«3 + Ti&n 2 + n 2 £ u O! + r 2 r 4 ^i 

+ f 4^1 + T 2 H^ + T 3 6^i + T$H 5 , (B.4) 

W v u cn = ^Nl + £ 2 ATf + F<h^i + JVifTgSi + F0i O2 S 2 + N 2 H& 2 . (B.5) 



After integrating out the heavy messenger fields we end up with the non-renormalizable 
operators as discussed in section [3l cf. also Figures HHSJ 

In addition to the renormalizable operators discussed so far there are six more operators 
allowed by the symmetries which are 

W™ = Tir 9 fii + T 2 r 9 3 + 3i3 3 & + T 9 ^ 2 + r 4 fi 2 ft 3 + rxHxH, . (B.6) 

The first two operators contribute effectively to the T\T 2 H^^ l operator already present and 
for the sake of simplicity we have not shown them in Figure EJ The third operator generates 
the dimension six operator FT 2 H^H 2 4<j) 2 (j) 23 which gives a contribution to the 2-2 element of 
the down-type quark and charged lepton Yukawa matrix. In fact the correction has the same 
phase and the same SU{5) Clebsch-Gordan coefficient as the leading order coefficient so that 
we can safely neglect it. The last three operators finally give, after integrating out the heavy 
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Figure 5: The supergraphs before integrating out the messengers for the up-type quark sector. 
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messengers, dimension seven and eight operators which give only small corrections (in our 
model we have discussed operators up to dimension six). The dimension seven operators, for 
instance, are induced by r^if^ which gives corrections to the 1 — 3 and 2 — 3 elements of 
the up-type quark Yukawa matrix which are very small compared to all other elements which 
are generated at maximum by a dimension five operator. 



References 

[1] J. Beringer et al. [Particle Data Group Collaboration], Phys. Rev. D 86 (2012) 010001. 

[2] P. Adamson et al. [MINOS Collaboration], Phys. Rev. Lett. 107 (2011) 181802 
[arXiv:1108.0015] . 

[3] Y. Abe et al. [DOUBLE-CHOOZ Collaboration], Phys. Rev. Lett. 108 (2012) 131801 
|arXiv:1112.6353| . 

[4] K. Abe et al. [T2K Collaboration]. larXiv:1106.28221 

[5] F. P. An et al. [DAYA-BAY Collaboration], Phys. Rev. Lett. 108 (2012) 171803 
|arXiv:1203.1669] : Y. Wang, talk at What is v! INVISIBLES 12 and Alexei Smirnov 
Fest (Galileo Galilei Institute for Theoretical Physics, Italy, 2012); available at 
|http://indico.cern.ch/conferenceTimeTable.py?confId=195985 

[6] J. K. Aim et al. [RENO Collaboration], Phys. Rev. Lett. 108 (2012) 191802 
[arXiv: 1204.0626] . 

[7] M. Ishitsuka, talk at Neutrino 2012 (Kyoto TERRS A, Japan, 2012); available at 
http ://kds.kek.jp/conferenceTimeTable.py?confId=9151 

[8] P. F. Harrison, D. H. Perkins and W. G. Scott, Phys. Lett. B 530 (2002) 167 
[hep-ph/0202074] . 

[9] H.-J. He, F.-R. Yin, larXiv:1104.2654l Z.-Z. Xing, larXiv:1106.3244l N. Qin and 
B. Q. Ma, Phys. Lett. B 702 (2011) 143 |arXiv:1106.3284] : Y. j. Zheng and B. Q. Ma, 
larXiv: 1106.40401 S. Zhou, larXiv: 1106.4808] T. Araki, larXiv:1106.52TTl N. Haba, 
R. Takahashi, larXiv: 1106. 59261 D. Meloni, larXiv:1107.022Tl S. Morisi, K. M. Patel 
and E. Peinado, larXiv: 1107.0696!; W. Chao, Y.-J. Zheng, larXiv: 1107.07381 H. Zhang, 
S. Zhou. larXiv: 1107. 10971 X. Chu, M. Dhen, T. Hambve. larXiv: 1107. 15891 P. S. B. Dev, 
R. N. Mohapatra, M. Severson. [arXiv:1107.2378l R. d. A. Toorop, F. Feruglio, C. Hage- 
dorn, larXiv: 1107.34861 Q. H. Cao, S. Khalil, E. Ma and H. Okada, larXiv: 1108.05701 
S. F. Ge, D. A. Dicus and W. W. Repko. larXTv:1108.0964l F. Bazzocchi. larXiv: 1108. 24971 

[10] S. Antusch and V. Maurer, Phys. Rev. D 84 (2011) 117301 [arXiv: 1107.37281 [hep-ph]]. 

[11] D. Marzocca, S. T. Petcov, A. Romanino and M. Spinrath, JHEP 1111 (2011) 009 
|arXiv:1108.06T4l [hep-ph]]. 

[12] S. F. King, C. Luhn and A. J. Stuart, Nucl. Phys. B 867 (2013) 203 |arXiv:1207.574Tl 
[hep-ph]]; C. Hagedorn, S. F. King and C. Luhn, Phys. Lett. B 717 (2012) 207 
|arXiv: 1205.3TT41 [hep-ph]]; I. K. Cooper, S. F. King and C. Luhn, JHEP 1206 (2012) 



22 



130 [arXiv: 1203.13241 [hep-ph]]; I. de Medeiros Varzielas and G. G. Ross. larXiv: 1203.66361 
[hep-ph] . 

[13] A. Meroni, S. T. Petcov and M. Spinrath, Phys. Rev. D 86 (2012) 113003 
[arXiv:1205.524T1 [hep-ph]] . 

[14] N. Haba, A. W atanabe and K. Yoshioka, Phys. Rev. Lett. 97 (200 6) 041601 
[hep-ph/0603Tl6] ; X. G. He and A. Zee, Phys. Lett. B 645 (2007) 427 [hep-ph/06071~63] ; 
W. Grimus and L. Lavoura, JHEP 0809 (2008) 106 [arXiv:0809.0226] : H. Ishi- 
mori, Y. Shimizu, M. Tanimoto and A. Watanabe, Phys. Rev. D 83 (2011) 033004 
[arXiv: 1010.3805] : Y. Shimizu, M. Tanimoto and A. Watanabe, larXiv:1105.2929l 
X. G. He and A. Zee, larXiv:1106.4359l 

[15] C. S. Lam, Phys. Rev. D 74 (2006) 113004 [hep-ph/06li017] ; C. H. Albright and 
W. Rodejohann, Eur. Phys. J. C 62 (2009) 599 [arXiv:0812.0436] : C. H. Albright, 
A. Dueck and W. Rodejohann, Eur. Phys. J. C 70 (2010) 1099 [arXiv: 1004.2798] . 

[16] S. F. King and C. Luhn, JHEP 1109 (2011) 042 |arXiv:1107.5332l [hep-ph]]: I. K. Cooper, 
S. F. King and C. Luhn, JHEP 1206 (2012) 130 [arXiv: 1203. 13241 [hen-ph]]: C. Hagedorn, 
S. F. King and C. Luhn, Phys. Lett. B 717 (2012) 207 |arXiv:1205.3TT4l [hep-ph]]. 

[17] S. Antusch, S. F. King, C. Luhn and M. Spinrath, Nucl. Phys. B 856 (2012) 328 
[arXiv: 1108 .42781 [hep-ph]]. 

[18] S. F. King, Phys. Lett. B 439 (1998) 350 |hep-ph/ 9806440] ; S. F. King, Nucl. Phys. B 562 
(1999) 57 [hep-ph/9904"21Q| ; S. F. King, Nucl. Phys. B 576 (2000) 85 [hep-ph/99 12492] ; 
S. F. King, JHEP 0209 (2002) 011 |hep-ph/0204360] ; S. Antusch, S. Boudjemaa, 
S. F. King, JHEP 1009 (2010) 096 [arXiv: 1003 .5498] . 

[19] S. F. King, JHEP 0508 (2005) 105 [hep-ph/0506297] . 

[20] S. F. King, Phys. Lett. B 675 (2009) 347 |arXiv:0903.3"T99] . 

[21] K. Harigaya, M. Ibe and T. T. Yanagida, Phys. Rev. D 86 (2012) 013002 |arXiv:1205.2l98l 
[hep-ph]]. 

[22] S. Antusch, S. F. King, C. Luhn and M. Spinrath, Nucl. Phys. B 850 (2011) 477 
[arXiv:1103.5930l [hep-ph]]. 

[23] A detailed review and references can be found in: G. C. Branco, R. G. Felipe and 
F. R. Joaquim, arXiv:1111.5332 [hep-ph]. 

[24] G. G. Ross, L. Velasco-Sevilla and O. Vives, Nucl. Phys. B 692 (2004) 50 
|hep-ph/0401064| . 

[25] S. F. King and M. Malinsky, Phys. Lett. B 645 (2007) 351 |hep-ph/0610250] . 

[26] F. Feruglio, C. Hagedorn and R. Ziegler, larXiv:1211.5560l [hep-ph]; M. Holthausen, 
M. Lindner and M. A. Schmidt. larXiv:1211.69~53l [hep-ph]. 

[27] H. Georgi and C. Jarlskog, Phys. Lett. B 86 (1979) 297. 



23 



S. Antusch and S. F. King, Phys. Lett. B 631 (2005) 42 [hep-ph/0508044] . 



S. Antusch, S. F. King and M. Malinsky, Nucl. Phys. B 820 (2009) 32 [arXiv:0810.3863l 
[hep-ph]]. 

G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo and A. M. Rotunno, Phys. 
Rev. D 86 (2012) 013012 |arXiv:1205.5254l [hep-ph]]. 

H. Minakata, H. Nunokawa, W. J. C. Teves and R. Zukanovich Funchal, Phys. Rev. D 
71 (2005) 013005 [hep-ph/0407326] ; A. Band yopadhyay, S. Ch oubey, S. Goswami and 
S. T. Petcov, Phys. Rev. D 72 (2005) 033013 |hep-ph/0410"283] . 



S. Antusch, S. F. King, M. Malinsky and M. Spinrath, Phys. Rev. D 81 (2010) 033008 
[arXiv:0910.5T27l [hep-ph]]. 

S. Antusch and M. Spinrath, Phys. Rev. D 79 (2009) 095004 |arXiv:0902.4644l [hep-ph]]. 
S. Antusch and M. Spinrath, Phys. Rev. D 78 (2008) 075020 [ arXiv: 0804 .0717 1 [hep-ph]]. 
A. Juttner [on behalf of the FLAG Gollaboration]. larXiv:1109.1388l [hep-ph]. 

G. Aad et al. [ATLAS Collaboration], Phys. Lett. B 716 (2012) 1 [arXiv:1207.72L4l 
[hep-ex]]; S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 716 (2012) 30 
[arXiv: 1207.72351 [hep-ex]]. 

U. Ellwanger, C. Hugonie and A. M. Teixeira, Phys. Rept. 496 (2010) 1 [arXiv:0910.1~785"1 
[hep-ph]]. 

J. C. Callaghan and S. F. King. larXiv:1210.6913 l [hep-ph]. 

S. Antusch, J. Kersten, M. Lindner, M. Ratz and M. A. Schmidt, JHEP 0503 (2005) 
024 [hep-ph/0501272] . 

K. Kowalska, S. Munir, L. Roszkowski, E. M. Sessolo, S. Trojanowski and Y. -L. S. Tsai, 
larXiv:1211.1693l [hep-ph]. 

Z.-z. Xing, H. Zhang, S. Zhou, Phys. Rev. D77 (2008) 113016. [ arXiv:0712.14l9l [hep- 
ph]]. 

K. Nakamura et al. (Particle Data Group), J. Phys. G 37 (2010) 075021. 

H. Leutwyler, Nucl. Phys. Proc. Suppl. 94 (2001) 108-115. | hep-ph/0011049] . 



24 



